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Abstract 



Realistic predictions are made for lepton flavour violating and flavour changing 
- - neutral current decays in the large tan/5 regime of a SUSY-seesaw model. Lepton 

^ ' flavour violating neutral Higgs decays are discussed within this framework along 

O ■ with the highly constraining charged lepton decays li —>■ ij'-f. In the b — s system 

I the important constraint from 6 — > 57 is studied in conjunction with the rare FCNC 

Bg fifi{Tfi) decays and the recently measured AM^. 
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C:? . 1 Introduction 
> 
o 

I Supersymmetric theories contain a number of possible sources of lepton flavour violation. 
\^ • As a result, the bounds on rates of LFV and FCNC are particularly restricting upon the 
O ■ SUSY parameter space. Even in the case of minimal supergravity(mSUGRA), where soft 
^ I SUSY breaking masses and trilinear couplings are flavour diagonal at the GUT scale, the 
presence of right-handed neutrinos and the PMNS mixings are enough to radiatively induce 
large LFV rates. It is these LFV rates, induced by renormalisation group (RG) running that 
we wish to study in the present work. 

D ! A large value of tan (3 enhances LFV and FCNC effects in SUSY theories and it is 

' interesting to study them together within a SUSY-GUT scenario. In the large tan/5 limit 
it has been noted that Higgs mediated contributions to FCNC decays may be dominant. 
I The most interesting example of Higgs mediated decays is that of —> /i+/i~[2, 1] in the 
;h ■ MSSM. This decay is proportional to tan^ /5 and so is greatly enhanced by large values of 
5^ ! tan (3. If the Higgs boson is very light then this rate could become very exciting. Another 
interesting case is the Higgs mediated contribution to the recently measured B^—B^ mixing 
at the Tevatron. Along with improved bounds on leptonic Bg decays these measurements 
add a complimentary b — s constraint to that of 6 — > S7. 

We propose to study the rates of lepton flavour violating decays of r, 5° and MSSM 
neutral Higgs bosons in a SUSY-seesaw model with bi-large neutrino mixing. From this 
framework we hope to make general predictions for some of the interesting lepton flavour 
violating decays mentioned above and also to study their correlations. The present work 
differs from previous studies in that it makes use of a top-down analysis of a Grand Unified 
SUSY model. Also we are studying lepton flavour violating Higgs decays in conjunction with 
Higgs mediated decays, such as Bg yU"''/^"- 

"'^Based on a talk given by J. Parry at the Mini-workshop on Particle Physics Phenomenology. Tainan, 
R.O.C., June 5th-6th 2006 and the article [1] 
^e-mail: jkparry@cycu.edu. tw 
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2 A large tan/? SUSY-seesaw model 



In the present work we shall study the MSSM+i/r constrained at the GUT scale by minimal 

SU{5) unification. In Minimal SU{5) the matter superfields of the MSSM+z/r are contained 
within the representations; 10 = {Q, If^, E'^), and 5 = (L, D'^) and 1 = {N'^). In addi- 
tion there are 5 = He) and 5 = (if^, He) Higgs representations. Using these matter 
superfields we can construct the Yukawa section of the superpotential as follows, 

K""^'^ = I io,io,h + 10,5,-H + y;^- 5,i,-h + ^ Ml la,-. (i) 

The first term of eq. (1) gives rise to a symmetric up-quark Yukawa coupling F„ and the 
second term provides both down-quark and charged lepton Yukawa couplings. As a result 
we have the GUT scale relation, Yf. = Yj. The final two terms of eq. (1) are responsible for 
the neutrino Yukawa coupling, Y^ and the right-handed neutrino Majorana mass, Mr which 
combine to produce the light neutrino mass matrix via the type I seesaw mechanism, 

mLL = -vlY,M^'Yj. (2) 



We begin by rotating the 10 and 5 representations such that Yd is diagonal. As a result 
Yu is solely responsible for the CKM mixing. The lepton sector is more complicated due to 
the structure of the see-saw mass matrix eq. (2). Hence, in this basis we have, 

Ya = Yf'^ 

Y, = UYf'^ (3) 
Mr = WM'^i'^W'^. 

Here the singlet neutrino has been rotated such that the combination, YjY^, is diagonal. In 
the case were, W — 1, the mixing matrix U is identified as the PMNS matrix^. Such a case 
may occur naturally, for example in models with a U{l)x fiavour symmetry. For simplicity 
we assumed W = 1 and that the above matrices are all real. 

Throughout the analysis it was assumed that both the neutrino Yukawa and Majorana 
matrices have a hierarchical form with Y^^^^ ~ 1 and Mr^ ~ 10^"^ GeV, as shown in eq. (4) 
and (5). We choose such a hierarchy in the neutrino Yukawa sector with one eye onb — t — T 
Yukawa unification and 5'O(10) with tan/3 ~ 50. 

y^diag < y^diag < y^diag ^ ^ (4) 

Mjf < M^;'^ < M^^'s ~ 10^^ GeV . (5) 

This hierarchy results in a neutrino mass ordering which is of the Normal Hierarchical type 
with, ^ 0.55 x 10"^ eV, m^^ ^ 0.009 eV, m^g ^ 0.05 eV. 

The consequence of our choice of X^^ ~ 1 is that it shall induce, via RG evolution, large 
fiavour mixing in the slepton mass squared matrix. This slepton fiavour mixing leads to 
a lepton-slepton mis-alignment and produces lepton fiavour violation in the charged lepton 
sector. The choice of a large neutrino yukawa coupling enhances this effect and will result 
in a large rate for r — > /I'j. It should be noted that such a large neutrino yukawa coupling 
is inevitable in models based on 5'O(10). 

In the soft SUSY breaking section of the Lagrangian fiavour blind boundary conditions 



^Such an SU{5) model has previously been used to study fj, — > 67, hadronic EDMs, etc. [4] 
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are assumed at the GUT scale. In this work we are interested in the renormahsation group 
effects of off-diagonal elements of the Yukawa couphngs as discussed above. Due to the 
large mixings in the PMNS matrix, these off-diagonal elements can be large even in the case 
of GUT scale flavour blindness. Throughout our analysis we have also assumed that the 
trilincar coupling, Aq = 0, at the GUT scale. Under these assumptions the soft Lagrangian 
becomes flavour blind and hence contains no new lepton flavour violating sources. 

3 LFV and FCNC Phenomenology 

Observations of neutrino oscillations imply the existence of massive neutrinos with large solar 
and atmospheric mixing angles. The small neutrino masses are most naturally explained via 
the seesaw mechanism with heavy singlet neutrinos. Even in a basis where both Yg and 
Mji are diagonal in flavour space Y^, is always left as a possible source of flavour violation. 
In SUSY models this flavour violation can be communicated to the slepton sector through 
renormahsation group running. The initial communication is from running between the GUT 
scale and the scale of Mji. Although the scale Mr is far above the electro- weak scale its 
effects leave a lasting impression on the mass squared matrices of the sleptons. Subsequently 
flavour violation can enter into the charged lepton sector through loop diagrams involving 
the sleptons and indeed such effects have been used to predict large branching ratios for 
T — > /X7 and —>■ e^j within the MSSM [5, 6, 7, 8]. As such, the experimental bound on 
charged lepton flavour violation provide important constraints on the possible flavour mixing 
of the MSSM and other extensions of the SM. 

In the Standard Model FCNCs arc absent at tree-level and only enter at 1-loop order. In 
extensions of the SM, the MSSM for example, there also exist additional sources of FCNC. 
A clear example comes from the mixings present in the squark sector of the MSSM. These 
mixings will also contribute to FCNCs at the 1-loop level and could even be larger than 
their SM counterparts. An example that wc shall study in this work is the flavour changing 
couplings of neutral Higgs bosons and the neutral Higgs penguin contribution to such decays 
as Bg fj,^lj,^ and B^ — B^ mixing. When examining such flavour changing in the quark 
sector it is important to consider the constraint imposed by the decay 6 — > 57. 

3.1 Higgs mediated FCNC 

It has been pointed out that Higgs mediated FCNC processes could be among the flrst signals 
of supcrsymmetry[2]. In the MSSM radiatively induced couplings between the up Higgs, 
and down-type quarks may result in flavour changing Higgs couplings. In turn this will lead 
to large FCNC decay rates for such process as Bg and 5° — i?^ mixing. In the 

Standard Model the predicted branching ratio for Bg — > A*^A*~ is of the order of 10~^, but 
in the MSSM such a decay is enhanced by large tan (3 and may reach far greater rates. It is 
also possible to extend this picture to the charged lepton sector where similar lepton flavour 
violating higgs couplings are radiatively induced. These LFV couplings can then result in 
such Higgs mediated decays as r — 3/i, the neutral Higgs decays , — > r/x and even 

combined with the afore mentioned FCNC couplings in the decay Bg ^ Tji. 

Loop diagrams induce flavour changing couplings of the kind, b'^sH^*. Similar diagrams 
with Higgs flelds replaced by their VEVs will also provide down quark mass corrections and 
will lead to sizeable corrections to the mass eigenvalues [9] and mixing matrices [10]. Fur- 
thermore the 3-point function and mass matrix can no-longer be simultaneously diagonalised 
[11] and hence beyond tree- level we shall have flavour changing Higgs couphngs in the mass 
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eigenstate basis. Such flavour changing Higgs couphngs can be summarized as, 



■'FCNC 



= -d 



Ri 
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RL 
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(6) 



These flavour changing couplings can in fact be related in a simple way to the finite non- 
logarithmic mass matrix corrections [3], 
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(7) 



where, Ago = (sq-/?, Cq_^, —i), for 5° = {H°, /lO, A^). It is clear that the FCNC couplings 
are related as, [Xrl] = [X^jK In general we should also notice that, [XjilI-j ^ [XRL]ji- 
Hence, in the case of, (i, j) = (6, s), we have [Xkl],,, ^ ^ [^lr^s- 



In the MSSM with large tan (3 the dominant contribution to Bg — > comes from the 
penguin diagram where the dilepton pair is produced from a virtual Higgs state [2]. Thus 
in combination with the standard tree-level term Cuh = ^VeJ-R.^LH'^ -\- h.c. the dominant 
tan P enhanced contribution to the branching ratio turns out to be. 
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where the matrix (5m^ is in the 1-loop mass eigenstate basis. 

In the large tan/5 limit the Higgs Double Penguin(DP) contribution to B^—B^ mixing is 
dominant [12]. Following the notation of eq. (6), we can write the neutral Higgs contribution 
to the AB = 2 effective Hamiltonian as. 



ayAB=2 
^eff 
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(9) 



where we have defined the operators. 
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It is common at this point to notice that the contribution to Q2 is dominant over 
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due to a suppression from the sum over Higgs fields, T 

contribution to Qi receives a factor JF+ = ( -^j^ + 
typically results in a suppression factor, ^ < T~ jT^ < Recalling that \Xi^r 
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-^[XjiLlbs, it may be possible for the Qf^^ contribution to give a significant effect. On the 
other hand, the contribution to Qf^^ is highly suppressed. 

Following the above conventions we can write the double penguin contribution as, 

- -^MbJeA 2^ Ml 

s s 

^ A/T -f2 r>LR RL]bs[X ];^ji]bs 

-^MnJsA 2. Ml 

s ^ 

Here Pf^^ — —1.06 and — 2.56, include NLO QCD renormahsation group factors 
[12]. After taking into account the relative values of JF^, the two P's and the factor of 2 
in eq. (13), we can see that there is a relative suppression, | < AM^^/AM^^ < |. This 
relative suppression shall be discussed further in the following section. 

There is a large non-perturbative uncertainty in the determination of /s^. Two recent 
lattice determinations provide [13, 14], 

^ 230 ±30 Me V (14) 

/2f = 259±32MeV, (15) 

which in turn give different direct Standard Model predictions for AMf^, 

^^SM'o4 ^ I7.8±8ps-i (16) 
^^SM'05 ^ I9.8±5.5ps-i (17) 

The recent precise Tcvatron measurement of AMg is consistent with these direct SM predic- 
tion but with a lower central value [?, 16] , 

AMf = 17.311^:1 ± 0.07ps-^ (18) 



3.2 Higgs mediated lepton flavour violation 

We have seen that flavour changing can appear in the couplings of the neutral Higgs bosons 
and is enhanced by large tan/5. In the quark sector, interactions of the form, dadiH^*, are 
generated at one-loop and at large tan f3 can become comparable to the tree-level interaction, 
dRd^H^. Similar Higgs-mediated flavour violation can also occur in the lepton sector of SUSY 
seesaw models through interactions of the form eRCLH^* . This leads to the possibility of 
large branching ratios for Higgs-mediated LFV processes such as Bg — > r//, r — > 3// and 
lepton flavour violating Higgs decays [17, 18, 19]. 

As we did for the down-quarks, we can find effective fiavour changing Higgs vertices in- 
duced after heavy sparticles are integrated out of the lagrangian. In the 1-loop corrected 
eigenstate basis these couplings are again related to the charged lepton finite mass correc- 
tions, (5me)^7*V^;„. 

Possibly the most interesting application of such lepton fiavour violating couplings is 
in the decays of MSSM Higgs bosons. These lepton fiavour violating Higgs couplings shall 
facilitate fiavour violating Higgs decays. We can write the partial widths of the lepton fiavour 
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violating Higgs boson decays within the MSSM as follows, 
1 {5m. 



167r 



- k^|'M5(i 



)^ {l-{x, + x,y){l-{x,-x,y) (19) 



= [sin(a-/?), cos(a-/3), -i] for S ^ [ , h° , ] . 



Here, — hP, represent the three physical Higgs states, Ms = Mh, M^, Ma are their 
masses, £j — t, fj,,e are the three charged lepton states and Xi — {m^jMsY- 

We can also make use of the LFV Higgs couplings to study the Higgs mediated contribu- 
tions to the process r ^ 3/i for example. The dominant Higgs contributions will come from 
the Higgs penguin contribution. In the MSSM with large tan the dominant contribution 
to the branching ratio of r — > 3// turns out to be. 



Br(T S/i) 



40 967r3 



Ml, 



(20) 



Here is the lifetime of the tau lepton and is the Yukawa coupling of the muon. The 
present experimental bound for this decay is as follows [20], 



BR(t ^ 3//) < 1.9 X 10-^ at 90% C.L. 



(21) 



Going one step further, the lepton flavour violating Higgs couplings can also be combined 
with the quark flavour changing coupling studied earlier. In this way we can also study the 
LFV and FCNC process Bg r/i. In the MSSM with large tan/3 the dominant Higgs 
contribution will again come from the penguin diagram mediated by Higgs bosons. The 
branching ratio for this decay may be written as. 
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(22) 



here Xi = {m^jM, 



Bs 



For the analysis two-loop RGEs for the dimensionless couplings and one-loop RGEs for 
the dimensionful couplings were used to run all couplings down to the scale M^r where the 
heaviest right-handed neutrino is decoupled from the RGEs. Similar steps were taken for the 
lighter M2R and Mm scales, and finally with all three right-handed neutrinos decoupled the 
solutions for the MSSM couphngs and spectra were computed at the Z scale. This includes 
full one loop SUSY threshold corrections to the fermion mass matrices and all Higgs masses 
while the sparticle masses are obtained at tree level. A detailed description of the numerical 
procedure can also be found in reference [1]. A function 'Y^i {Xf^ — X'^^'^Y / is evaluated 
based on the agreement between the theoretical predictions and 24 experimental observables. 

The analysis uses approximately 100 evenly spaced points in the SUSY parameter space 
each with values of itlq and M1/2 in the plane, itlq = 300 — 1400 GcV and M1/2 = 300 — 1000 
GeV. For each of the 100 points we also held fixed /i = +120 GeV, tan/3 = 50 and Aq = 0. 
The remaining input parameters are allowed to vary in order to find the minimum of our 
function. 
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4 Discussion 



In this section we shall present the results of our analysis. We shall discuss the LFV and 
FCNC phenomenology is turn. 
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Figure 1: A contour plot of Br(r fj,-f) x 10^ in the mo — M1/2 plane with fi = 120 GeV 
and ^0 = 0. The present experimental bound is marked by the contour labelled "3.1". 



Let us firstly look at the important lepton flavour violating decay ii ij'y. A contour 
plot showing the variation of the branching ratio for t fi'y in the mo — M1/2 plane is given 
in fig. 1. This plot shows that the rate for this process is very high over most of the plane 
and even in the low mo — M1/2 region the rate can exceed the present experimental upper 
bound. These large rates are a direct consequence of our choice of Y^^^^ ~ 1. 
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Figure 2: The Higgs mediated contribution to the branching ratio for Bg ^^{t^j) 
is plotted in the lower (upper) panel against the Pseudoscalar Higgs mass, m^o. The 
horizontal line at 0.8 x 10"'' shows the present experimental upper bound. 



The plot for the Higgs mediated contribution to the flavour changing neutral current 
decay, Bg — fi~^iJ,~, is presented in the lower panel of fig. 2. As previously reported, the 
branching ratio for Bg — >■ is particularly interesting with rates ranging from 10~^° for 
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heavy m^o up to almost 10^^ for a particularly light rriAO- The present 90% CL. experimental 
upper bound of 0.8 x 10~^[21] excludes the highest predicted rates and hence is beginning 
to probe the Higgs sector into the region m^^o ~ 150 GeV. This rare decay is certainly of 
particular interest and future studies will continue to probe the Higgs sector to an even 
greater extent. Within the standard model the expectation is that Bt{Bs — > fJ^'^fJ^~)sM ~ 
10^^. Therefore we can see from fig. 2 that the Higgs mediated contribution to this process 
would dominate as long as m^o < 500 GeV and this decay has the possibility of being the 
first indirect signal of supersymmetry. 



Figure 3: The correlation of Bt{Bs and AM^^ are given for fs^ = 230 MeV 

(upper panel) and /b^ = 259 MeV(lower panel). The horizontal line in both panels shows 
the present experimental bound on Br(i?s — > The vertical lines mark the central 

value of (AMf - AMf*^). 



The lepton flavour violating decay Bg Tjj, plotted against the Pseudoscalar Higgs 
mass, is shown in the upper panel of fig. 2. This plot shows that the branching ratio for 
this decay could be as large as ~ 10~^°, but are certainly not as interesting as the flavour 
conserving decay just discussed. The branching ratio for the decay r — >■ 3/i are of a similar 
order of magnitude to those in the upper panel of fig. 2. 

In the limit of large tan/3, and AMg are correlated. This correlation is 

shown in the two panels of fig. 3. For these two panels the two different values of Jb^ listed 
in eq. (14,15) are used. The upper panel (/^^ = 230 MeV) shows that the central value 
of (AM^^^ — AMf^) coincides with the bound from Bt{Bs —>■ ix^ix~). The lower panel 
{fss = 259 MeV) shows that the data points with AMf ^ at the central value, are in fact 
ruled out by the bound on Bt{Bs — >■ The uncertainty in the SM prediction for 

AMg is rather large and in fact all of the data points of fig. 3 are allowed by the recent 
tevatron measurement. These two panels clearly show that the interpretation of the recent 
measurement depends crucially on the uncertainty in the determination of Jb^- 

The plot in fig. 4 shows the ratio, AM^^ /AM^^', of the contributions to the operators 
Qf^^ and (^2^ as defined in eq. (13). It is commonly assumed that the contribution to the 
(^2^ operator, AM^^, is dominant. From fig. 4 we can see that the contribution of the Qf^^ 
operator, AM^^, is between 40% and 90% of AM^^ and hence is significant. 

Our numerical results for the branching ratio of r/z are presented in fig. 5c. The 
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rates for this process are plotted against the Pseudoscalar higgs mass, tjiao- The predicted 
rates range from 10"^ up to a few times 10~^. In fig. 5c we can see a broad band of points 
stretching from m^o = 100 GeV up to m^o = 600 GeV showing that most points arc 
predicting a rate of around 10"^. The decay rate appears to be almost independent of the 
Pseudoscalar higgs mass, although there is a slight peak at the lower range of the higgs mass, 
which is where the highest rates are achieved. Notice that the data points shown in fig. 5 are 
divided into two groups. The grouping is determined by whether the Bg fi^fi~ bound is 
in excess or not, as indicated in the figure caption. From this grouping we can see that the 
points with the highest predicted — > t/j, rates appear to be excluded by the Bs — * 
bound. In this way the Higgs mediated contribution to the Bs decay can provide additional 
information on the allowed Higgs decay rate. 
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Figure 4: Plot for the ratio, AM^^ /AMg^, of the Higgs contribution to the operators 
Qf^^ and Q^^ 



Fig. 5b shows our predictions for the lepton flavour violating Pseudoscalar higgs decay, 
^4° — > TfjL. The rates for the decay of the Pseudoscalar are almost identical to those of the 
heavy CP-even higgs shown in fig. 5c. 

The rate for the decay of the lightest higgs boson is shown in fig. 5a where it is again 
plotted against the Pseudoscalar higgs mass. The predicted rates for this decay show a 
very different dependence upon m^io, as they are spread over a large range from 10^^^ to 
10~^. The branching ratio for the hghtest higgs appears to be inversely proportional to the 
Pseudoscalar higgs mass. Hence this LFV decay will only be interesting if m^o < 300 GeV 
where its rate can be comparable to those for the other neutral higgs states. The data points 
for the lightest Higgs boson decay of fig. 5a have again been divided into two groups. As 
before, the two groupings depend upon whether the Bg — >■ /^"""/^^ bound is being exceeded 
or not. The plot shows that the data points for particularly light Pseudoscalar Higgs mass 
are excluded by the Bg — > A'-"''/^" bound. These excluded points correspond to the largest 
predictions for the decay /i° — > t/jl and leaves, 3 x 10~^, the highest allow decay rate. 



5 Conclusions 

We have analysed a supcrsymmetric-seesaw model constrained by SU (5) unification at the 
GUT scale. We have been concerned with making predictions for lepton flavour violating 
decay processes such as r — ^'-f, ^ ^ e'j and H^, Tfi. With lepton flavour violation 

in mind, we chose to study the large tan (3 region of parameter space where such effects are 
enhanced. In addition we chose Y^^^^ ~ 1 under the influence of SO{10) unification in which 
tan /3 is naturally large. 
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Our numerical procedure utilises a complete top-down global fit to 24 low-energy 
observables. Through this electroweak fit we were able to analyse lepton flavour violating 
rates for rare charged lepton decays. The choice of a large third family neutrino Yukawa 
coupling naturally leads to large 23 lepton flavour violation. Within this scenario we have 
also used these numerical fits to analyse the FCNC processes Bg fifi^Tfi), 5° mixing 
and the lepton flavour violating decays of MSSM Higgs bosons. As a result, we have been 
about to make realistic predictions for this phenomenology while ensuring ii ij'-f and 
b ^ sj constraints are satisfied. 

Our model predicts a r — >■ /i7 rate near the present experimental limit. We have seen 
that the branching ratio for — > r/i could be interesting with a rate as high as few 10~^. 
The Higgs mediated contribution to the branching ratio for Bg fi^fi^ was found to be 
particularly large and could even exceed the current experimental bound. In addition to the 
constraint of r — >■ n'y, we also saw that the Bg — >■ bound may also act as a stringent 

restriction on the allowed rates for 0° Tfj,. We found the rates for the LFV decays Bg 
and r ^ 3yU are rather small < 10^^*^. 
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Figure 5: BR(ff'^ Tjj.) plotted against the Pseudoscalar Higgs mass, m^o. For each 
of these fitted points we have assumed the following values, = 120 GeV and Aq = 0. 
Here the dots(circles) mark points for which Br{Bg — > is below(above) the present 

experimental limit. 



The correlation of Bg —>■ fj^^fi'^ and AMg in the large tan/3 limit was also studied. The 
constraint from the recent tevatron measurement is highly dependant upon the determination 
of J'bs ■ It was also found that the Higgs contribution to the operator Qf^^ can compete with 
that of the (^2^ operator. 
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